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INTRODUCTION 


The  first  thing  one  learns  about  sound  in  the  ocean  is  that  it  fluctuates.  Such  varia- 
bility is  much  more  than  merely  conspicuous;  it  is  quite  general  and  ever  present.  Thus, 
variability  exists  in  all  parameters  of  the  sonar  equation  such  as  propagation  loss,  target 
strength,  ambient  noise  and  recognition  differential,  though  each  time  due  to  rather  differ- 
ent causes.  Fluctuations  in  the  ocean  are  also  relevant  to  understanding  the  performance 
of  naval  acoustic  equipment  in  that  they  unmistakably,  and  sometimes  drastically,  influence* 
such  performance.  Typical  is  the  observation  that  targets  tend  to  be  detected  at  ultra -long 
ranges  during  periods  of  signal  surges  and  lost  at  short  ranges  during  periods  of  signal 
fades . To  put  it  differently,  the  actual  detection  probability  falls  off  more  slowly  in  field 
exercises  than  would  be  expected  in  the  absence  of  fluctuations. 

Of  the  many  kinds  of  variability  mentioned  above,  fluctuations  in  propagation  loss  have 
received  a great  deal  of  attention.  In  consequence,  a substantial  amount  of  information  both 
experimental  and  theoretical  has  been  gathered  concerning  candidate  mechanisms  of  trans- 
mission fluctuability.  It  has  thus  become  apparent  that,  under  certain  conditions  at  least, 
sound  wave  scattering  off  the  temperature  microstructure  of  the  medium  can  account  for 
most  of  the  observed  variability.  Hereafter,  we  shall  refer  to  this  mechanism  and  the., 
attendant  body  of  theory  as  the  Kolmogorov,  Chernov,  and  Tatarski  (KCT)  framework/ 

In  this  paper  we  offer  a straightforward  generalization  of  the  KCT  results  away  from 
the  transversal  correlation  restriction  inherent  therein.  The  academic  value  of  this  work 
is  limited;  the  tactical  value  thereof  is,  however,  paramount  in  that  for  real  life  employ- 
ment of  sonar  detection  the  transversal  correlation  case  is  of  vanishing  measure. 

The  observation  fundamental  to  our  work  is  that  the  ocean  is  usually  in  a state  of  tur- 
bulent motion.  Correspondingly,  the  values  of  the  temperature  at  every  point  in  the  ocean 
undergo  irregular  fluctuations;  the  values  of  the  temperature  at  different  spatial  point  at 
the  same  instant  in  time  also  differ  from  one  another  in  a random  fashion.  What  has  been 
said  applies  as  well  to  all  other  oceanographic  quantities.  In  particular,  since  the  index  of 
refraction  of  the  ocean  is  a function  of  temperature  and  salinity,  we  can  take  the  viewpoint 
that  the  refractive  index  field  is  random.  To  become  specific,  we  shall  assume  that  the 
Kolmogorov  theory  of  locally  homogeneous  and  isotropic  turbulence^)  provides  a sufficiently 
good  description  of  the  refractive  index  microstructure . 

To  extract  information  concerning  the  randomness  of  an  acoustic  wave  propagating 
through  this  turbulent  and  unbounded  ocean,  we  make  use  of  the  wave  equation  to  connect 
the  statistical  properties  of  the  random  medium  to  the  implied  statistical  properties  of 
wave  parameters  within  the  framework  of  a correlation  theory.  We  accomplish  this  only 
to  first  order  in  perturbation  theory  thus  restricting  the  realm  of  validity  of  our  results  to 
high  frequencies  and  small  refractive  index  fluctuations. 
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The  main  output  of  this  analysis  is  the  structure  functions  of  the  logorithmic  amplitude 
and  phase  fluctuations.  For  conciseness,  we  have  explicitly  handled  only  the  amplitude 
fluctuations,  but  the  same  treatment  can  be  continued  with  no  difficulty  to  include  the  phase 
fluctuations  as  well. 

We  thus  find  that  the  correlation  length  of  the  logorithmic  amplitude  fluctuations  is 
given  by, 
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where  9 measures  the  angle  between  the  observation  plane  and  the  wave  front,  X is  the 
wave  length,  L the  distance  between  the  source  and  the  receiver,  Lg  the  outer  scale  of 
turbulence,  and  0^  is  the  solution  of 

/xr  = 4l,Lo 

cos  9 X sin  0 
We  also  find  that, 

<^(log  A/A0)2])=  .31  C2  k7/6  L11/6 
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index  fluctuations,  k stands  for  the  wave  number  2tt/x  and  where  A/Aq  represents  the 
amplitude  of  the  perturbed  wave  normalized  to  the  plane  wave  amplitude  Ag  . 

The  first  three  chapters  are  intended  to  familiarize  the  reader  with  the  mathematical 
language  of  random  fields  extensively  used  herein,  with  the  physical  picture  of  turbulent 
flow,  and  with  the  corresponding  microstructure  of  the  refractive  index  field  respectively. 
The  wave  equation  is  developed  in  chapter  D.l,  and  the  small  wavelength  modification  of 
Rytov  to  the  first  order  perturbation  theory  appears  in  D.2.  We  solve  the  Rytov  equation 


using  the  method  of  spectral  expansions  In  chapter  D.3.  Explicit  evaluation  of  the  solution 
for  the  Kolmogorov  model  is  included  In  chapter  D.  4.  The  analysis  concludes  with  a chap- 
ter on  comparison  with  experimental  data (3)  gathered  under  conditions  where  the  turbulence 
mechanism  is  expected  to  dominate . 

A.  RANDOM  FUNCTIONS 

• A random  function  on  T is  a family  of  random  variables  {§(t),  ?(s), . . . - 
corresponding  to  all  elements  (t,  s,  . . .)  in  the  set  T . 

• We  shall  regard  the  random  function  |(t ) as  specified  if  for  each 

subset  (t.,  t„,  t ) in  T we  are  given  the  distribution  function 

pt1,t2>...,tn  (W-’V  sp{5(t1><x1.---.;(tn)<yn}  . (A-l) 


In  most  actual  physical  problems  it  is  excessively  complicated  to  measure  these  distribu- 
tion functions . They  are  also  too  cumbersome  to  be  used  in  practice.  It  is,  therefore, 
convenient  to  restrict  oneself  to  the  simplest  numerical  characteristics  of  the  multidimen- 
sional distribution,  the  moments 


m,  m. 


Hiyiiy...  2<^  ^1^  ? (t2^  * 


m 

• ? <V> 


(A-2) 


Here, 


M,1  “to  m 

<?  (tp  ? ct2) ...  ? n(tn)>  s 


(A-3) 


r mi  m2  m 
S J X1  X2  •••XnndF  , (X,,X9,...,x  ) 

-a>  11  1 u2  * * Ln  -L  ^ n7 


where  the  integration  on  the  right-hand  side  is  to  be  performed  in  the  Stieltjes  sense.  If 
the  distribution  F has  a density, 
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the  moment  reduces  to  the  standard  form: 
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We  shall  further  restrict  ourselves  by  taking  into  account  only  those  properties  of  a random 
function  which  are  determined  by  its  first  and  second  moments,  the  mean  value  and  the  corre- 
lation function: 


m(t)  = <|(tr)  > ; B(t,s)  = <5(t)  ?(s)  > . <A"6> 

• The  random  function  % (t)  is  called  stationary  if  all  the  finite-dimensional 

distribution  functions  defining  % (t)  are  invariant  under  translations  of  the  whole  group 

of  points  (t, , t_, . . . , t ) in  T , 

12  n 


Fti+T  t2+T,...,tn+T  <X1,X2.....Xn)  Ft15t2, 

If  |(t)  is  stationary,  then  obviously: 

m(t)  = m ; B(t,s)  = B(t-s) 
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The  physical  content  of  the  concept  of  stationarity  is  clear.  It  means  that  ) describes 
the  time  variation  of  a numerical  characteristic  5 of  an  event  such  that  none  of  the  ob- 
served macroscopic  factors  influencing  the  occurrence  of  the  event  changes  in  time. 


• To  determine  the  mean  value  m(t)  and  the  correlation  function  B(t,  s)  of  a 
random  function  £(t)  we  must  first  take  a large  number  N of  realizations  of  E;(t),  written 

and  then  calculate  the  arithmetic  mean  of  ?(j)  (t)  over  j 
for  every  value  of  t , or  the  arithmetic  mean  of  ^ ^ ^(t)^^(s ) for  every  pair' of  values 
t and  s . However,  in  practice,  observation  of  a random  function  and  the  subsequent 
processing  of  the  data  usually  turn  out  to  be  quite  complicated.  It  would,  therefore,  be  very 
desirable  to  be  able  to  get  along  with  as  small  a number  of  realization  as  possible.  Indeed, 
the  practical  value  of  the  correlation  theory  of  stationary  random  functions  is  to  a consid- 
erable extent  due  to  the  fact  that  if  £(t ) is  stationary,  its  mean  value  m and  its  corre- 
lation function  B(t)  can  usually  be  calculated  by  using  just  one  realization  of  |(t)  . This 
is  the  content  of  the  ergodic  theorem. 


More  precisely,  if  5 (t ) is  a stationary  random  function  satisfying  certain  quite  general 
conditions,  then 
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where  the  integrals  are  defined  as  the  limits  of  the  corresponding  approximating  sums: 
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We  notice,  that  in  the  statement  of  the  ergodic  theorem,  the  limit  is  taken  to  mean  limit  in 
the  mean -square. 

• A stationary  random  function  can  be  represented  in  the  form  of  a stochastic 
Fourier-Stieltjes  integral, 


CO 
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(A-13) 


where  stationarity  requires: 
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Hence : 


00 


B(t-s)  = j 


(A -15) 
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becomes, 


00 


B(t-s)  = J du,elu)(t>SVa) ) . 


(A -16) 


Obviously,  W(od)  ^0  , and  therefore  the  Fourier  transform  of  a correlation  function  must 
be  non-negative.  Khinchin  showed(^)  that  the  converse  assertion  is  also  true:  if  the  Fourier 
transform  of  B(t)  is  non-negative,  then  there  exists  a stationary  random  function  ^(t) 
with  B(t)  as  its  correlation  function.  W(<»)  is  called  the  spectral  density  function  of  §(t ) . 


• In  general,  the  random  functions  encountered  in  practical  applications  can  very 


often  be  regarded  to  a high  degree  of  accuracy  as  stationary.  However,  the  opposite  can 
also  occur.  The  mean  value  of  some  oceanic  or  meteorological  variable  undergoes  com- 
paratively slow  and  smooth  changes,  and  hence  they  are  random  functions  that  cannot  be 
regarded  as  stationary. 

In  the  case  where  ?(t)  represents  a non -stationary  random  function,  i.e.,  where 
< § (t)>  changes  in  the  course  of  time,  we  can  consider  instead  of  |(t)  the  difference, 


For  values  of  t which  are  not  too  large,  slow  changes  in  the  function  |(t ) do  not  affect 
the  value  of  this  difference,  and  it  can  be  regarded,  at  least  approximately,  as  a stationary 
random  function: 


Random  functions  that  possess  this  property  are  called  random  functions  with  stationary 
increments  (Kolmogorov) . (2)  It  is  easy  to  show  that  the  correlation  function 


AT(t)  = ?(t+T)  - 5(t)  . 
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< A (t )>  =-  0 . 

T 


(A -18) 


is  expressible  as  a linear  combination  of  the  structure  functions: 

D(t,s)  =<(?(t)  - f-(s))  > . (A -20) 

For  B^(t,s ) to  depend  only  on  the  difference  (t-e)  , it  shall  suffice  to  have, 

D(T)  = < (§(t+T)  - 5(t ))2>  . (A -21) 

Roughly  speaking,  the  value  of  D(  t)  characterizes  the  intensity  of  those  fluctuations  of 
% ( t ) with  periods  which  are  smaller  than  or  comparable  with  t . 

• A random  function  with  stationary  increments  can  be  represented  in  the  form, 


?(t)  = ?(0)  + J (1  - ela)t)d0U) 


(A -22) 


where  |(0)  is  a random  variable,  and  the  Fourier -Stieltjes  amplitudes  d0(u))  obey  the 
condition 

<d0( u))d0*(  \)>  = 6( or  X)W(u))dujdX  . (A-23) 


Hence, 


D(t)  = 2 J d«)  (1  - cos  wt)W(a>) 


(A -24) 


(A -25) 

B.  THE  MICROSTRUCTURE  OF  TURBULENT  FLOW 

For  us,  the  most  important  fact  about  the  ocean  is  that  it  is  usually  in  a state  of  tur- 
bulent motion.  We  shall  therefore  need  some  basic  information  concerning  the  statistical 
properties  of  developed  turbulent  flow.  The  statistical  theory  of  turbulence,  was  initiated 
in  the  papers  of  Friedmann  & Keller.  A very  important  advance  was  achieved  in  1941, 

when  Komogorov,  ^ Obukhov^  and  later  Onsager,  von  Weisacker  and  Heisenberg 


and 


W(«)  ^ 
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established  the  laws  which  characterize  the  basic  properties  of  the  microstructure  of 
turbulent  flow  at  very  large  Reynolds  numbers. 


Consider  an  initially  laminar  flow  of  a viscous  fluid.  This  flow  can  be  characterized 
by  the  values  of  the  kinematic  viscosity  v , the  characteristic  velocity  scale  v and  the 
characteristic  length  L . The  quantity  L specifies  the  dimensions  of  the  flow  as  a whole 
and  arises  from  the  boundary  conditions  of  the  fluid  dynamics  problem.  Suppose  now,  that 
for  some  reason  or  other,  a velocity  fluctuation  occurs  in  a region  of  size  A of  the 
basic  laminar  flow.  The  characteristic  time  r^=  A/v^  which  corresponds  to  this  fluctua- 


tion specifies  the  order  of  magnitude  of  the  time  required  for  the  occurrence  of  the  fluctua- 
tion. The  energy  per  unit  mass  of  the  given  fluctuation  is  given  by  v^2  . Thus,  when  the 
velocity  fluctuation  under  consideration  occurs,  the  amount  of  energy  per  unit  time  per 
unit  mass  which  goes  over  from  the  initial  flow  to  the  fluctuational  motion  is  equal  in  order 
of  magnitude  to, 


i 


Vk3/A  • 


On  the  other  hand,  the  local  velocity  gradients  thus  developed  are  given  by  the  ratio  v' ^ A 
and  therefore  the  energy  dissipated  as  heat  per  unit  mass  of  the  fluid  per  unit  time  is  of 
order  of  magnitude 


C = v vJ^2/A2 


If  the  velocity  fluctuation  which  arises  is  to  have  existence  of  its  own,  it  is  clearly  neces- 
sary that  the  inequality 

i 

I 

v’a3/A  > \|V^2/a2  (B-l) 

hold,  i.e.,  that 

Re  A = A v’/  v > 1 . (B-2) 

Since  all  these  calculations  are  accurate  only  to  within  undetermined  numerical  factors,  it 
would  be  more  appropriate  to  write  instead  of  equation  B-2. 

Re , > Re 
A cr 

Where  Re  ^ denotes  the  "inner"  Reynolds  number  corresponding  to  fluctuations  of  size  A , 

and  Re  is  some  fixed  number  that  cannot  be  determined  precisely.  Notice,  however, 
that  if 


Re  . > Re 
A cr 


> 


the  direct  energy  dissipation  for  the  velocity  fluctuation  with  size  A is  small  compared  to 
the  energy  it  receives  and  thus  the  fluctuation  can  transfer  almost  all  its  energy  to  smaller 

perturbations.  Consequently,  the  quantity  v '^ / A , which  represents  the  energy  per  unit 


mass  received  per  unit  time  by  eddies  of  nth  order  from  eddies  of  (n-l)'th  order  and 
thereupon  transferred  to  eddies  of  (n+l)th  order,  is  a constant  for  perturbations  of  almost 
all  sizes.  In  the  smallest  velocity  perturbation  with  size  /u  , this  energy  is  converted 

U 2 2 

into  heat.  The  rate  of  dissipation  of  energy  into  heat  is  of  the  order  € ~ vVg  / Ag  > 
and  hence  for  velocity  fluctuations  of  all  scales  we  have 


€ ~V^/A 


(B-3) 


for  the  energy  transmitted  down  the  eddies  ordered  chain.  It  is  naturally  clear  that  if  the 

condition  vL/\i  > Re  is  not  met  for  the  flow  as  a whole,  the  laminar  motion  is  stable. 

cr 


The  largest  eddies  which  arise  as  a result  of  the  instability  of  the  basic  flow  are  of 
course  not  isotropic,  since  they  are  influenced  by  the  special  geometric  properties  of  the 
flow.  However,  these  special  properties  no  longer  influence  the  eddies  of  sufficiently  high 
order,  and  therefore  are  good  grounds  for  considering  the  latter  to  be  isotropic.  Then, 
since  eddies  with  dimensions  much  larger  than  | p | do  not  influence  the  two  point  function 
jj7  (r  + "p)  - v (r ) J'  this  difference  will  depend  only  on  isotropic  eddies. 


We  thus  arrive  naturally  at  the  scheme  of  a locally  isotropic  random  field.  The  random 
field  is  hence  uniquely  characterized  by  the  longitudal  structure  function  D ( p)  . The  form 
of  the  structure  function  can  be  established  by  using  the  qualitative  consideration  developed 
above.  In  fact,  consider  a value  of  p that  is  large  compared  to  the  inner  scale  Ag  of  the 
turbulence  and  small  compared  with  the  outer  scale  L-o  of  turbulence.  Then  the  velocity  dif- 
ference at  r and  r + p is  mainly  due  to  eddies  with  dimension  comparable  to  p . But  the 
only  parameter  which  characterizes  such  eddies  is  the  energy  dissipation  rate  * . Thus,  we 
can  assert  that  D ( p)  is  only  a function  of  p and  $ . Dimensional  arguments  then  provide 


D/r  ( p)  = C(cp)2/3 


(B-4) 


Ag  « p « Lq 

where  C is  a dimensionless  constant  of  order  unity.  This  equation  represents  the  famous 
Kolmogorov  "two-thirds  law."  For  p « Ag  , changes  of  velocity  occur  smoothly,  since 

now  the  relative  motions  are  laminar.  The  velocity  difference  can  therefore  be  expanded 
in  power  series  of  p , and 
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D ( p)  = a p2 

it 


p«A0 


(B-5) 


FIG.  1.  GENERAL  SHAPE  OF  THE  STRUCTURE  FUNCTION  D ( o ) 

ti 

In  equation  (B-4)  it  is  assumed  that  p « Lq  , where  Lq  is  the  outer  scale  of  the 
turbulence.  When  p is  increased,  the  condition  p « Lq  is  violated.  Then  the  large 
eddies,  which  cannot  be  regarded  as  isotropic  and  homogeneous,  begin  to  influence  the 
value  of  v(r+  p)  - v(r ) . In  this  case,  the  structure  function  D ( p)  depends  on  the 

n 

coordinates  of  both  observation  points,  and  no  universal  law  can  be  given  which  describes 
the  structure  function  for  large  values  of  p . We  can  only  state  that  the  growth  of  the 
structure  function  slows  down  for  p » Lq  . Figure  1 shows  the  general  shape  of  the 

longitudinal  structure  function  D ( p)  . For  small  values  of  p , the  curve  can  be  re- 
placed by  a parabola  with  great  accuracy,  then  the  part  of  the  curve  corresponding  to  the 
"two-thirds"  law  begins,  and  finally,  in  the  region  of  the  outer  scale  of  turbulence,  the 
curve  starts  to  saturate. 
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In  many  applications,  it  is  often  expedient  to  regard  the  values  of  p as  not  being 
bounded  from  above  by  the  value  of  Lq  . However,  if  we  cannot  neglect  the  saturation  of 

the  structure  function,  it  is  necessary  to  use  interpolation  formulas  which  approximately 
describe  the  behavior  of  the  structure  function  for  large  values  of  p . For  small  values 
of  p , these  formulas  must  reduce  to  the  same  values  of  D (p)  as  given  in  equation  (B-5). 
In  what  follows  we  shall  use  the  Karman  function, 


D (p) 

// 


!-2 


■\}  - hxm  (j£)  ki/3  (■%)  ] • 


(B-6) 


,.2 


Here  <v'  > is  the  mean  square  velocity  fluctuation,  Lq  the  outer  scale  of  turbulence  and 

^1/3  (li^  ) *s  t^ie  order  Bessel  function  of  the  second  kind  of  imaginary  argument. 

For  -=£-  « 1 the  Karman  structure  function  is  approximately  equal  to 
L0 

2/3 


Vrt-sp&sj  <v'2>(t£) 


(B-7) 


and  hence  coincides  for  p « Lq  with  the  Kolmogorov  "two- thirds"  law. 

C.  TURBULENT  MIXING  OF  PASSIVE  ADDITIVES 

• Although  the  foregoing  discussion  shall  prove  quite  essential,  let  us  observe 
that  it  is  really  the  structure  of  the  temperature  and  salinity  fields  that  ultimately  deter- 
mine the  characteristics  of  wave  propagation  in  the  free  ocean,  in  that  the  wave  velocity  C 
is  a function  of  temperature,  salinity  and  depth.  We  shall  now  take  the  view  that  the  tem- 
perature and  salinity  fields  can  be  regarded  as  conservative,  passive  additives  to  the  tur- 
bulent ocean.  If  a volume  of  liquid  is  characterized  by  a concentration  9 of  additive,  then 
by  saying  that  the  additive  is  conservative  we  mean  that  the  quantity  9 does  not  change  when 
the  volume  element  is  shifted  about  in  space.  By  the  additive  being  passitive  is  meant  that 
the  quantity  9 does  not  affect  the  dynamical  regime  of  the  underlying  turbulence. 

As  usual,  we  separate  the  value  of  the  quantity  9 into  the  mean  value  < 9 > and  the 
fluctuation  0'  , 

9 = <B>.+  9’ 

and  define  the  following  measure  of  the  inhomogeniety  of  the  spatial  distribution  of  9 , 
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(C-l) 


G = | Jdv<0,2> 

V 

so  that  G vanishes  if,  and  only  if,  0’  vanishes  identically  in  V . Using  the  equation  of 
molecular  diffusion  which  must  be  obeyed  by  the  concentration  0 of  the  passive  additive 
we  can  easily  show  that, 

||  = J dV  [k(v<0  >)2  - D <(v  0'  )2  >]  (C-2) 

V 

where  D represents  the  coefficient  of  molecular  diffusion  and  K the  coefficient  of  turbu- 
lent diffusion,  with  K»D  . For  clarity,  we  recall  that  if  represents  the  mean  flow  of 
0 caused  by  molecular  diffusion, 

qM  = - dv  < 0 > 

and,  similarly,  if  represents  the  density  of  turbulent  flow  of  0 , 

Op  = - Kv<  0 > . 

Thus,  the  time  dependence  of  G is  governed  by  the  delicate  interplay  of  two  distinct  phys- 
ical processes,  turbulent  mixing  and  molecular  diffusion.  As  a result  of  turbulent  mixing 
the  inhomogeneity  of  the  spacial  distribution  of  8 is  increased  and  large  local  gradients 
of  0 are  created.  Only  after  these  large  gradients  have  appeared  does  the  process  of 
molecular  diffusion  play  a significant  role  by  smoothing  out  the  spatial  distribution  of  0 . 

The  quantity, 

N = D<(v  0'  )2>  (C-3) 

represents  the  amount  of  inhomogeneity  which  disappears  per  unit  time  due  to  molecular 
diffusion  and  it  is  analogous  to  the  energy  dissipation  rate  e . Correspondingly,  the 

2 

quantity  K(v<  0 >)  represents  the  amount  of  inhomogeneity  which  appears  per  unit  time 

due  to  the  turbulence,  and  is  similar  to  v A > the  rate  of  production  of  the  energy  of 
the  velocity  fluctuations . 

We  can  carry  out  an  even  more  detailed  analogy  between  the  velocity  fluctuation  in  a 
turbulent  flow  and  the  concentration  fluctuations  of  a passive  additive  0 . Concentration 
inhomogeneities  0’  with  geometrical  dimensions  A appear  as  a result  of  the  action  of 
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velocity  field  perturbations  with  dimensions  A and  characteristic  velocities  v A . The 

amount  of  inhomogeneities  appearing  per  unit  time  due  to  turbulence  is  clearly  given  by, 

<2 

VA  ® a / ^ • The  rate  of  levelling  out  of  the  inhomogeneities  0*  A is  of  the  order 

of  D0*  ,2/  A2  . When 
A 

vA  0'  A2/  A » D0’  A2/ A2  (04) 

A A A 


that  is,  when 

AvA/D  » 1 , (05) 

the  inhomogeneity  0*  A is  not  dissipated  by  the  action  of  molecular  diffusion,  but  rather 

has  a stable  existence  and  can  subsequently  subdivide  into  smaller  eddies.  This  process 
of  subdivision  proceeds  until  inhomogeneities  appear  for  which 

Ava/D  ~ 1 . (06) 

These  inhomogeneities  are  dissipated  by  the  process  of  molecular  diffusion  at  a rate  equal 
to  N . 

Thus,  the  amount  of  inhomogeneity  transferred  per  unit  time  from  the  largest  eddies 
down  the  chain  to  the  smallest  eddies  is  constant  and  equal  to  the  rate  N at  which  the  in- 
homogeneity is  dissipated.  Again,  the  largest  inhomogeneities  in  the  distribution  of  0 are 
not  isotropic.  However,  the  smallest  inhomogeneities  can  be  considered  isotropic.  Hence, 
0(r+p)  - 0(r ) , determined  mainly  by  inhomogeneities  of  size  p , can  be  considered 
statistically  isotropic  for  values  of  p « Lq  . It  follows,  then,  that  0(r)  can  be  regarded 
as  a locally  isotropic  random  field.  Correspondingly,  the  structure  function  D(  p)  can 
depend  only  on  p,  N,  e . Dimensional  considerations  lead  to. 


D(  P)  = C02 


(C-7) 


Aq  « p « h) 


which  is  the,  by  now  familiar,  "two-thirds"  law  for  the  concentration  of  a passive  additive 
(Obukhov).  (6) 


D.  PARAMETER  FLUCTUATIONS  OF  ACOUSTIC  WAVES  PROPAGATING 
IN  A TURBULENT  OCEAN  AT  SMALL  WAVELENGTHS 

1 • The  Wave  Equation  for  Acoustic  Propagation 

Consider  a fluid  medium  in  a state  of  equilibrium.  We  shall  denote  the  equilibrium 
density  by  and  the  equilibrium  pressure  by  PQ  . The  actual  density  at  the  point  r 
and  at  time  t shall  be  denoted  by  p(r>t)  . Hence,  the  relative  change  in  density  due 
to  the  passing  of  the  acoustic  wave  will  be  given  by 

p (r,t)  = p(r,t)  - pQ  . 


Similarly  the  difference  between  the  actual  pressure  and  the  equilibrium  value,  will  be 
denoted  by, 

A 

p(r,t)  = p(r,t)  - pQ  . 


We  must  now  find  an  equation  giving  the  dependence  of  p on  r and  t . To  accomplish 
this,  we  observe  that  the  total  pressure  p and  density  p satisfy, 

dv 

p = -vp  (Newton's  second  law)  (D-l) 


+ V(  pv)  = 0 


(continuity  equation) 


(D-2) 


where  _v  represents  the  velocity  of  the  acoustic  oscillations.  Discarding  terms  of  second 
and  higher  orders  of  smallness  in  the  small  quantities  p,  p,  and  v , one  has: 


„ _ - 
% at  " 


(D-3) 


+ pd  = o 


(D-4) 


Differentiating  with  respect  to  time  the  continuity  equation  D-4 


2 * 


(D-5) 
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We  can  now  replace  dv/dt  in  equation  D-5  from  Newton's  equation  D-3  to  obtain: 


2 a 

o o „2  * 

— 5 - 7 p 

at 


= o 


(D-6) 


A A 

To  find  one  more  relationship  between  p and  p , we  must  make  an  assumption  about  the 
thermodynamic  character  of  the  acoustic  process.  Specifically,  we  shall  assume  that  alter- 
ations in  pressure  and  density  are  so  rapid  that  no  heat  energy  has  time  to  flow  away  from 
the  compressed  part  of  the  fluid  before  this  part  is  no  longer  compressed.  Such  compressions 
are  said  to  be  adiabatic.  Evidently,  the  amount  of  heat  changes  by, 


(D-7) 


and  introducing  the  specific  heats 


we  have: 

dQ  ■ cp($ )p  dv  + <=»  (f  )v  dp  • 


(D-8) 


If  the  fluid  is  a perfect  gas,  pV  = RT  , and  hence, 
dT  _ dp  , dV 

T ■ T — • 

providing: 


With  these  relations,  the  heat  change  equation  becomes: 


(D-9) 


(D-10) 
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For  adiabatic  processes,  dQ  = 0 and  hence  equation  D-10  provides, 
dp  _ dV 

T " “ Y~v 

where 

= Cp 
Y " CV 

We  rewrite  equation  D-ll  in  the  form, 


(D-ll) 


(D-12) 


P = 


or,  rather, 


P = 


for, 


C = 


(D-13) 


If  the  fluid  is  a liquid,  the  relationship  between  p and  p given  in  equation  D-13  remains 
valid  but  the  constant  C2  no  longer  relates  to  equilibrium  values  in  exactly  the  same  way. 
Differentiating  twice  with  respect  to  time,  equation  D-13  becomes: 


Eliminating  with  the  aid  of  equation  D-6,  the  density  p 


, 2 * 0 
^ ^ - v2p  = 0 
<T  at: 


(D-14) 


which  is  the  desired  wave  equation.  The  constant  C , that  might  in  principle  depend  on 
r , represent  the  wave  velocity. 
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2.  The  Rytov  Method  of  Approximation 


To  describe  the  propagation  of  an  acoustic  wave  through  a turbulent  ocean,  we  must 
solve  the  wave  equation 


v2p  - 


2 

•2-f=o 

at 


_4 

for  an  index  of  refraction  N (r ) that  is  a random  function  of  r . Notice,  that  here  we 
have  regarded  the  distribution  of  inhomogeneities  as  static,  neglecting  their  change  as  a 
result  of  heat  conduction,  diffusion,  convection  and  drift.  We  can  neglect  this  change  if 
the  propagation  time  is  small  compared  to  the  characteristic  time  scale  of  change  in  the 
inhomogeneities.  Correspondingly,  if  the  source  is  monocromatic,  the  time  dependence 

of  p(r,t ) is  of  the  typical  form  elu)  C and  the  wave  equation  reduces  to  the  well  known 
Helmholtz  equation. 


V2p  + k2N2(r)p  = 0 


(D-15) 


where  k s ^ represents  the  wave  number  of  the  source, 
L0 


k 


We  set 


N(r)  = 1 + n(r ) 


with 


|n(r)  |« 


1 


and  initially  attempt  to  solve  the  Helmholtz  equation  by  the  method  of  small  perturbations. 
The  idea  behind  this  method  is  that  due  to  the  smallness  of  n(r)  the  solution  p(r)  differs 
little  from  the  non-turbulent  solution  and  hence,  that  p(r ) can  be  looked  for  in  the  form  of 
a power  series  expansion  in  some  small  parameter  of  order  n . One  can  think  of  the  mth 
term  in  the  expansion  as  representing  the  non-turbulent  wave  scattered  m times  by  the  tur- 
bulent inhomogeneities  of  the  medium,  and  hence  of  the  wave  actually  present  at  the  observa- 
tion point  as  the  sum  of  such  multiply  scattered  waves.  Indeed,  consider  the  non-turbulent 
wave  Pq(t  ) satisfying  the  homogeneous  Helmholtz  equation 


v2p0(r)  + k2p0(?)  = 0 
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and  allow  the  medium  to  develop  a turbulent  inhomogeneity  n in  the  index  of  refraction  at 
some  arbitrary  point  r Q . In  the  presence  of  Pq(p)  the  inhomogeneity  at  becomes 

the  source  of  a secondary  wave  p-^(r)  with  source  strength  -2k2n(r)Pg(r)  . The  sec 

ondary  wave  coexists  with  Pq(? ) and  satisfies  the  nonhomogeneous  Helmholtz  equation, 

v2p1(r)  + k2p1(r)  = - 2k2n(r)pQ(r)  . 


Consequently, 

pl(?)  = J d?0  e(?P0)  [-2k2n(rQ)p0(r0)] 


where  the  Green's  function  G(r  |r^)  satisfies 

V2G(r|rQ)  + k2G(r|r0)  = - 6(r-rQ) 

and  represents  the  contribution  to  the  wave  at  r due  to  a point  source  of  unit 
strength  at  . The  actual  field  measured  at  r is  therefore  given  by 

P(r)  = P0(r)  + J drQ  G(r|rQ)  r-2k2n(rQ )P0(rQ)]  • 

Let  now  another  inhomogeneity  develop  at  . In  the  presence  of  p-^(r)  , it  will  become 
the  source  of  a secondary  wave  of  strength  -2k2n(r  )p^(r  ) and  hence, 

P2(r)  = J drx  G(r  jr^)  ^-2k2n(r1)p1(r1)] 

By  repeating  the  argument  ad  infinitum,  it  is  not  difficult  to  show  that  the  wave  actually 
measured  at  r can  be  written  as  a multiple  scattering  expansion: 


p(r)  = pQ(r)  + £ Pm(r) 

m=l 


where 


[- 


2kn(r  . )p  _(r  . 

v m-1  *m-l  m-J- 
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£ The  iterative  nature  of  this  equation  indicates  that  the  mth  term  in  the  expansion  is  of  order 

nm  . Therefore,  for  sufficiently  small  n one  might  hope  that  higher  terms  in  the  multiple 
scattering  expansion  can  be  neglected  with  impunity.  Since,  however,  one  expects  the  fluc- 
tuations away  from  the  nonturbulent  value  to  grow  with  the  range,  a straightforward  applica- 
tion of  the  method  of  small  perturbations  shall  not  suffice  for  all  but  very  short  ranges. 

• Rather,  (?» 8)  we  write  equation  D-15  as, 

+ kV  = 0 

P 


and  notice  that  it  is  equivalent  to, 

v2  log  p + (v  log  p)2  + k2N2  = 0 . (D-16) 

We  set 

t 5 log  p 

so  that 

= log  A + iS 

where  A represents  the  amplitude  of  p,  S the  phase  thereof,  and  have, 

v2i)f  + (Vi|r)2  + k2^  = 0 . (D-17) 

We  are  now  ready  to  apply  the  method  of  small  perturbations  to  this  modified  differential 
equation  by  taking,  to  first  order, 

♦ = ^0  + *1  (D-18) 

where  the  zeroth  approximation 

*0  = log  s0  + is0 

satisfies: 

v%0  + (n0)2  + k2  = 0 (d-19> 
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and  where  ^ is  the  once  scattered  wave.  Then,  using  equations  D-l 7,  18,  19 

+ V^C  2v*q  + 7^)  + 2k  2n  + k2n2  = 0 (D"20> 


where 


*1  = log  ^ + i(S~V  * 

2 2 

In  equation  D-20,  we  omit  k n which  is  of  second  order  of  smallness  and  further  assume 
that 

W*il«H0l  • (D-21) 

Hence,  the  equation  to  be  solved  becomes: 

72^  + 2 v tJtq  ’ 7^  + 2k2n  = 0 . (D-22) 


Since, 

|v*0l~k 

condition  D-21  describing  the  limits  of  validity  of  our  approximation  now  reads: 

X|7ilr1|«  2tt  (D-23) 

and  expresses  the  smallness  of  the  change  of  ^ over  distances  of  the  order  of  a wave 

length.  This  inequality  implies  the  smallness  of  the  relative  amplitude  change  over  a 
wave  length 

\|v  log  tl\«  2tt  (D-24) 

A0 

and  the  smallness  of  the  change  of  the  phase  fluctuation  over  X , 

X|v(S-S0)|«  2tt  . (D-25) 
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The  amplitude  condition  D-24  Is  always  fulfilled  in  a weakly  inhomogeneous  medium, 
|n  |«  1 . The  phase  condition  D-25  implies  that  the  angle  of  inclination  of  the  ray 
the  initial  direction  is  small.  In  fact,  if 


we  have, 


for 


and  hence,  using  equation  D-25, 


Since  large  scale  inhomogeneities,  X - < Aq  , produce  sharply  directed  forward  scatter- 
ing, the  phase  condition  D-25  can  be  met  if  we  assume: 


3.  The  Method  of  Spectral  Expansions 

We  now  consider  the  problem  of  fluctuations  of  a monocromatic  wave,  confining  our- 
selves to  the  case  where  the  wave  length  X is  small  compared  to  the  inner  scale  of  tur- 
bulence Aq  . In  this  case,  the  angle  of  scattering  of  the  wave  by  refractive  index 
inhomogeneities  is  of  order  no  greater  than  X/Aq  and  is  thus  small.  Therefore,  the  value 
of  ^(r  ) can  only  be  appreciably  affected  by  the  inhomogeneities  included  in  a cone  with 
vertex  at  the  observation  point,  with  axis  directed  towards  the  wave  source  and  with  angular 
aperture  X/  A q . Hence  inhomogeneities  lying  on  the  opposite  side  of  the  source  from  the 
observation  point  never  contribute  to  scattering  and  we  shall  therefore  take  n(r ) = 0 in 
that  region.  To  simplify  the  calculation  we  shall  furthermore  assume  that  the  wave  is  plane. 
We  then  locate  the  origin  of  coordinates  on  the  boundary  of  the  region  occupied  by  the  refrac- 
tive index  inhomogeneities  and  direct  the  x-axis  along  the  direction  of  propagation  of  the 
incident  wave . 
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FIG.  2.  THE  DISTRIBUTION  OF  REFRACTIVE  INDEX  INHOMOGENEITIES 

Then, 

= (ik,  0,  0) 

and  the  differential  equation  D-22  reduces  to, 

9 S'h  2 

7 ^ + 2ik  — ^ f 2k  n = 0 . 


If  \«  Ag  > one  can  show  that  it  shall  suffice  to  solve  the  simpler  equation, 

2 d’h  o 

7 Ati  + 2i  + 2k'  n = 0 

where 


The  method  of  spectral  expansion  consists  in  taking 


(D-26) 


n(?)  = J e lx'r  dZ(x) 


(D-27) 
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and  searching  for  the  solution  in  the  form, 


♦■!_( ?)  = J eix'^  dcp  (*)  . 


(D-28) 


Using  the  differential  equation  D-26, 


i L ~ 2kx//')  dcp(n)+  2k^  V dZ(x)  - 0 


and  therefore, 


d®(H)  = ^ 


x /2k  + x - ic 

i.  // 


dZ  (it)  . 


(D-29) 


Notice  the  small  imaginary  part  that  has  been  added  to  the  denominator  to  render  the  expres- 

sion  meaningful  at  x = - h /2k  . Before  proceeding  any  further,  we  must  first  explicitly 

//  1 

incorporate  the  boundary  condition 
n(r ) =0  , X < 0 

into  the  Fourier-Stieltjes  transform  D-27  of  the  refractive  index  field.  We  therefore  take, 
n(r  ) = 0(X  ) Tj  (r) 

with  the  step  function  0(X)  defined  by 


0(X)  h 


1 , X > 0 

0 , X < 0 


and  invert  the  Fourier-Stieltjes  transformation  D-27  to  have, 


dZ(Z>  = - dK>  Jdv(co,  ?A)  J dXei('"-"i  >x 


(D-30) 
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where  d\j(x)  is  the  Fourier  Stieltjes  transform  of  Tl(r)  . 

Using  equation  D-30  we  can  now  formally  express  the  solution  / in  terms  of  the 
known  Fourier-Stieltjes  transform  of  Tl(r)  . Thus,  we  can  express  the  transform  of  the 
random  field  of  logorithmic  amplitude  fluctuations, 


Ref^r)  = J elK‘r  da(x) 

in  terms  of  dv(x)  using  the  obvious  relation 

da(x)  = ^ TdcoCn)  + dco*(-x)] 


and  obtain 


da(x)  = - tz  d*  |dv(cD,K,)  f 77 I 

4tT  " J X L xV2k  + H - ic  n 


dX  e 


i(«)-K  )X 


+ 


1 

/ 2k  - k + ic 

L // 


J dXe1^"*.*] 

0 


(D-31) 


We  can  also  express  the  Fourier-Stieltjes  transform  of  the  random  field  of  phase  # 

fluctuations 


. — ♦ — ♦ 

Im^Cr)  = J elH*r  d0(x) 

in  terms  of  d\>(x)  using, 

do(n)  = ^dcp(x)  - d<p*(-x)J  . 
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We  get, 


do(ir)=:  - rad\  Jd',<")>\>  [ 


x /2k  + H 

X n 


- 1 € 


- JdXe1*®-*.* 


0 


K /2k  - X + if  n 
X //  u 


J dXei(a,'\)x] 


(D-32) 


Correspondingly,  we  can  now  calculate  the  correlation  function  of  logorithmic  amplitude 
fluctuations  in  terms  of  the  known  correlation  function  of  the  refractive  index  field,  or  the 
Fourier  transform  thereof.  Indeed, 


Bft(p)  s < Reijf1(r+p)Re*iJr1(r)  > = 

= r e-iir,?  < da(it)da*(H'  )> 


and  recognizing  that 


< da(x)da*(x')>  = 6(x-x')  «(x  ,x?;x  )dx  dn'di?  dn' 

1 A n n X « u X i. 


we  have, 


B,(  p)  = [ dx  dx'dx  e 1 1 e 

R J //  n X 


ix  p in  (r  + p ) -ix'r 


' " e * ' Vk,’h1;hx)  • 0>33> 


The  spectal  density  function  n>K\',K  x)  *s  easily  connected  to  the  spectral  density 
function  of  the  refractive  index  field  via  its  definition  and  equation  D-31, 


?/2k+H  -i€)(H?/2k+HT+i€) 

X //  X n 


J 


i(cu-x)X  -i(® -x’  )X’ 
dXdX’e  " e " 
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(*?/ 2k+n  -i€)(H;/2k-H,-ic)  n 
La  1 * u 


r i(a)-H  )X  -iCtu-H^X' 

J dXdX'e  u e 


+ —5 r dXdX’e 

(Hf/2k-H  +i€)(H,/2k+K’+i€)  n 

J.  n La  U 


i( (o-K  )X  -i( (u-h*  )X' 


(D-34) 


_i f 

i -rVw2/9k-v  r'l  *' 


dXdX’e 


i(®-K  )X  -i( (u-k'  )X* 


(k  /2k-n  +i c ) ( h / 2k- h ’ - i c ) 


] 


where, 


<dv(to  ,k  )dv*(w  ' >h^)>  - 6( (u-u)' ) 6 (h^-k^)  $n  ( ® jH^dojd®  ’dn^d*^ 


We  are  now  in  a position  to  evaluate  the  h ,k’  integrals  in  equation  D-33.  To  do  so, 

an 

define : 

_ 7 ion  -iBH1 

F,  (a,  P;h  ) — dH  d*'  e "e  " $.(k  ,k’;k,)  . 

J-  . J n n A tt  n l 


Replacing  $ from  equation  D-34,  interchanging  the  (u)  ,X,X,')  integrals  with  the 

i\ 

) integration,  and  employing  repeatedly  the  spectral  representation  of  the  step 
function, 


SCO  = lim  |d, 

f -»0 

v m 


,-iojt 

u)+i  € 
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we  obtain: 


a B 

J dX  J dX’  Fn(X-X’;Hx)  * 


(D-35) 


0 0 


♦ 


where, 


00 


Fn(x-X';HX)  = J da,  e 


ioj(X-X' ) 


$n  ( 0)  > H ) 


— CD 


In  equation  D-35,  the  upper  limit  of  both  integrals  can  be  set  equal  to 
Y = Min  (a,  B)  • 

In  fact,  because  of  the  directional  characters  of  the  scattering,  the  waves  scattered  by  the 
layer  bounded  by  the  planes  X = Min  ( a > P)  and  X = Max  ( a » B)  are  incident  on  the  re- 
ceiver at  Max  ( a > B-)  but  not  on  the  receiver  at  Min  (a»B)  • Therefore,  these  waves 
can  be  neglected  in  the  calculation  of  the  correlation  function. 

Introducing, 


I = X-X’ 


T)  h X+X’ 


FA(a  ,B  ; hx)  becomes: 
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FA(a,P;nx)  = \ **  J dEF^?;^)  (y-  ?) 


H2  k2 

fees  (“ ^ ( B - a-  ?))  + cos  ( B - a+  f ) j ; 


|k2J  d?Fn(?;Hx)^  ^ 

n H 


(D-36) 


H2  K2 

[Sin  (~2k  ( B+  a'  o)  - sin  (“2^  ( B+  a-2  Y+  ?))]  . 

In  general,  the  function  Fn(  ? ; k ) falls  off  very  rapidly  to  zero  for  k £ 1 
fore,  the  important  contribution  to  the  values  of  the  integrals  occurs  for  g ^ — 

2 Hx 

1 . X.?  K 


There  - 
. In  the 


region 


r <; — 


, we  have 


much  less  than  the  inner  scale  of  turbulence  Aq  . But  Aq 


We  assumed  above  that  the  wavelength  X is 
■7 — , where  k _ is  the 


it  "x 

e of  l 

largest  wave  number  for  which  F ( ? ; h±)  still  differs  from  zero.  Therefore  we  have: 


±m 


im 


H 


jum 


x <_i r 


« 1 


H2? 


Thus,  « 1 in  the  important  region  of  integration  on  5 . 


Furthermore,  we  shall  be  interested  in  the  correlation  function  of  Rety  (r)  only  for 
values  of  the  argument  which  are  small  compared  to  the  range  y . This  means  that  in 

equation  D-36  we  consider  only  values  of  k which  satisfy  the  condition  — < v . Since 

1 x Hx 

% £ — in  the  important  region  of  integration,  within  this  region  we  have  5 « y . 
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Taking  all  these  simplifications  into  account,  we  obtain: 


x2  x2 

P^(o.>P>xJL)^k  !”y  cos  ( B~  o.)  J + —*>  sin  (~ 5^  ( 8 + O-  - 2y)) 

U Kx 

x2  v 

" 4 sin  (iK  (3+ a))]  Jd5F  (5;k  ) 

Hx  0 

Recognizing  that, 

V oo 

J d?Fn  (?;Hjl)  to  j d?Fn  (C;xx)  = n $n  (0;^) 

0 0 

we  finally  have: 


Now,  in  accordance  with  equation  D-33,  the  correlation  function  of  the  logarithmic  ampli- 
tude field  is  given  by, 


yj)  * j e 
«— ♦ 

Choosing  one  of  the  two  h 
have: 


lx  p 

x x Fa(r  +p  ,r  ;x  ) . 

axis  along  *p^  , we  can  perform  the  angular  integration  to 
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CO 


Vp)  a 2"  I ’VKVMP  WWV 


0 

We  are  therefore  interested  in, 

a h r +p  ; 3 = * 

//  //  " 

and  with  the  notation  of  figure  3 we  thus  have 
3-  a = p sin  0 ; 3 + a = 2L 


Px  = p cos  0 
V = L - p sin  9 


FIG.  3.  THE  GEOMETRY  OF  THE  CORRELATION  FUNCTION  MEASUREMENT 
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Then, 


Notice  that  B^(  p)  does  depend  on  the  range  L , thus  indicating  that  the  field  Re  ^(r  ) 

is  not  homogeneous.  This  circumstance  reflects  the  nonhomogeneous  choice  of  "boundaries" 
implicit  in  the  assumption  of  forward  scattering  and  is  not  characteristic  of  the  field  itself. 

9 4.  PHENOMENOLOGICAL  MODEL 

We  are  now  ready  to  choose  some  reasonable  model  for  the  Fourier  transform  $n(x.) 

of  the  correlation  function  characterizing  the  refractive  index  inhomogeneities.  A natural 
candidate  for  the  structure  function  of  a turbulent  field  has  already  been  suggested  in  chap- 
• ter  B.  Here,  however,  with  the  field  assumed  homogeneous  and  isotropic  rather  than 

locally  so,  the  structure  function  can  be  eliminated  in  favor  of  the  correlation  function 

B(  p)  = \ [d(»)  - D(  p)]  . 
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Using  the  Karman  proposal,  we  therefore  have, 

2/3  2 . Nl/3  . v 

Bn(  p)  " I W3)  <n  > (if;)  Kl/3  (l£) 


and  correspondingly, 


- r(ii/6) 

VK)  = ~375T 


n ' r(l/3) 


L°<n2>^^76 


In  what  follows,  w-e  shall  find  it  convenient  to  write 


VK>  " ^ 3/2^/  L0  <n2>  n4S2 l2)11/6 

y-CT  n r(l/3)  (1+k  Lq) 


(D-39) 


and  postpone  the  taking  of  the  limit  to  the  end  of  our  calculation.  With  this  choice  for 

* ( v ) » the  correlation  function  for  the  logarithmic  amplitude  becomes: 

nv 

ba(p)  = <n2>k3L^  [sVVA  PsineAi  + A3] 


(D-40) 


where 


H 

I *idHxJo(Hip  cos  0)  COS  ("2E  p sin  0 ) 7“'"272"ar/6" 
0 <1+KXV 


-y  k. 


A,  s 


3 J KxdKxJ0^Kxp  COS  ^ 
0 


sin  (■£  l) 


-YH, 


K. 


(1+^Lo)11/6 


A3  H J cos  0) 


sin  (-§£-  p sin  0 


) 


2 

e'YXi 
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We  begin  the  evaluation  of  recognizing  that 


A±  = Re  I(z) 


for 


I(z) 


■I 


H ' KxdKie 


-zh,  J.CSk,) 

X U 1 


,,  , 2.2.11/6 
(1+Kj.V 


and  where, 


z = y-i  a ; Rez  > 0 


a = p sin  0 


p s p cos  e 


Changing  variables,  we  can  write, 


I(z)  = -^2  J dx 


“zx/L0  Jo(^  x%) 


2I* 


(1+x) 


U75 


and  using  the  power  series  expansion 


Jo  (jl  x%)-  y tif{4)n*n 

KL0  ^=0  (no2 
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we  have: 


2 , „2.n 


I(Z)  = ?0  ^7?  © l dx 


-zx/l: 


,0(l+x)-^6  xn 


0 0 


(D-41) 


vn  , „2*n 


I(2)  ““2  I ^""FTT-  (^)  Y (n+1'  n+^  I 4 ) 

4L0  L0 


2L0  n=0 


where  Y (a  ( c jx)  represents  the  confluent  hypergeometric  function  of  Tricomi^ 

satisfying  the  integral  representation, 


Y(a|c|x)  = f dt  e"xt  ta_1  (1+t) 

T(a)  J 


c-a-1 


Rea  > 0 


Rex  > 0 


If  |JL|  <<  1 , we  can  use  the  behavior  near  the  origin  for  the  Tricomi  function, 


f ai|l4) 

L0 


r(5/6) 

r(ii/6) 


r(-5/6) 


i|r(  n+l  | n+1/6  | z/L2) 


r(n-5/6)  /z^/6_n 

rCn+1)  s) 


n ^ 


1 
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to  obtain, 


I(z) 


T(5/6)  1 

fnx/e)  ^2 


r(-5/6) 


+ 


(-l)n 

(n!)2 


F(n-  5/6) 


Employing  the  standard  integral  representation  of  the  Euler  gamma  function 


r( n-  5/6 ) = j dte_t  tn-11^6  , (n-5/6)  >0 

0 


I(z)  i niZiL.J^rc-B/G)  (4) 
rcn/6)  2I2  V 


5/6 


,5/6 


,n  , „2^sn 


2Itt  Lo 


(n!) 
n=l  v ' 


and  recognizing  that 


I (-#)"  ■ * [<-4f]  - 1 

n=l  ^n‘; 


« 
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it  follows 


1(2) 


1 H 5/6) 


r(' 


5/6)  (-4) 
L0 


5/6 


+ 2? 

2L0  L0 


5/6 


J 


dte-t  t’11/6 


(Jo[2(- 


M]  -o 


Consider  now, 


»(-£)  * J 


-t  . -ii/6  /,  vn(_£t 


2.  N^- 


dte  t 


KC2(-F)  ] -0 


and  use. 


J0[2(-t|T]  -1  “ ^ I dSJl(?C%) 


<4f 


to  have: 


<4f  . 

d(4)  . . j « J *.-*  t-e/6  „£(_£)*  ] 

0 0 


Notice,  however,  that 


dte-C  t'8'*  ™ e‘?2/4  | * (ll/6|21?2/4) 
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where  $(alc|x) 


represent  the  confluent  hypergeometric 


function  of  Kummer. 


(9) 


Hence, 


D 


( = T(l/6) 

v 4 z)  frn 


32/4z 

J d*  e'x 
0 


$ (ll/6|2|x) 


and  upon  using 


® [e'X  * (iVelUx)]  e"x  * (ll/6|2|x) 

we  find, 

2 2 

D(~fe)  = r<‘5/6)  [e"P  /4"  $ (ll/6|l|e2/4z)  -l]  . 


Thus: 


I(Z)  "*^72  f(  11/6)  + 775  r("5/6)(^)  e~B  ^ * (H/6|1|P2/4z) 
2I10  2L0  L0 


or, 


Hjl  «i 

l: 


where  use  has  been  made  of  the  Kummer  transformation, 
$(ajc|x)  = ex  f (c-a|c|-x)  . 
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The  real  part  of  equation  D-42  is  then  the  value  of  the  integral  we  have  been  trying 
to  evaluate . 

Consider  next, 


2 

-Y* 

Hx^Hxe 


sin(aHx  ) 


J0(e*x) 


Notice  that  G(  y , a ) has  been  so  chosen  as  to  provide 

A2  = G(y,  L/K)  ; A3  = G(y,  o/2k  sin  0) 

thus  reducing  both  A 2 andA^  to  an  evaluation  of  G(y}  a).  With 


. , 2.  1 

sin(  ax, ) 


— ^ = ol  j d Xcos  (Xan^) 

H.  n 


we  have: 


G(v?  a)  “ 


„l  "r  -(v-iox)x^  J0(BHp 

f M %d*xe 


0 0 


2 T 2 \H/ 6 
(i+X^  Lq) 


1 

+ f J dx  J H^e 
0 0 


-(v+iaX)xx  Jo(Pkx} 


2 T 2.11/ 6 

(i+xx  L0) 
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and  hence; 


1 

G(v,a)=f  J d\  TlCy-iaX)  + I(Y+iaX)  ] 

0 

Y+i  a 

G(v  ,a)  = ■£[  J dz  I(z) 
y-i  a 

G ( y > a ) = Im  K.  ( z ) 

where, 

K(z)  = J dz  I(z)  . 

Using  the  expression  for  I(z)  obtained  in  equation  D-42 

9 5/6 

K(z)  "*  \ Sj  (^)  + 7^5  n-5/6)  (^Ar) 

L0  2L0  4L0 

J dz  (-i§ )5/6  e~^/4z  $ (11/6  |l  j 02/4z) 

3 

Denoting, 

D(z ) = J dz  (-^|)  e~  ® ^4z  $ (ll/6|l|e2/4z) 

P 


it  follows  that. 


*«)  = - 4 I <4)  (■ 4)1V 6 ^ » ai/6|i|e2/4z) 


and  using 


[X"11/6  e’X  * (17/6 11  lx)] 


r(-5/6)  .-17/6  -x . % 

r(-ll/6)  X e * (11/611  lx) 


we  obtain: 


“ -4  » (-I)"11'6  * (17/6|1,B2/4z)  . 


Hence, 


(D-43) 


14  l «i 


Collecting  the  various  partial  results  developed  above,  we  find: 
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V •> 


2tt^  <r»2> 


,.3t3  r(ll/6) 

k L0  (1/3) 


lim  * 

y-0+ 


|^Im  K(y+i  ^ sin  9)  - sin  9 Re  ICy+i-^  sin  0) 


-Im  K(y+i^)  + ^ Re  ICy+i-^  sin 


4n  LQ  . 

p <<:  X sin  9 ’ <<:  L0 

For  normalization  purposes,  let  us  consider: 

Ba(0)  - 2tT*  <n  > k3L3  ^lim+  I(y)-  Im  K(y+i 


Upon  taking  3 = 0 in  both  I(z)  and  K (z)  we  obtain: 


I(V>"  HT176)  + ^5  n-5/6)  (-^) 


5/6 


irf  +i  L1  - 1 H5/6)  v+iL/k  1 , . f v+i  L/k\ 

x (y+i  ^ -(-n/6)  I n-11/6)  ^-5 ) 

*t)  L0 


11/6 


having  employed 

$ (a |c  |0)  = 1 . 
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Therefore, 


Ba(0-)  - < n2>  k3!3  I™  25'6  r(-ll/S)  sin  (^)  * 

(D-44) 

k-11/6  (a!*)-11/6  L11/6 
and  introducing  the  normalized  correlation  function 

bA(  p)  = ba(p)/ba(^) 
we  shall  have: 


for 


6 = (2Lq.) 


11/6 


[2 5/6  IX -11/6)  sin  tt)] 
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The  result  of  equation  D-44  is  quite  important  in  itself,  because 

where  o ^ stands  for  the  variance  of  the  logorithmic  amplitude  distribution.  In  fact,  we 
find: 

2 ^ 2^k7/6  t11/6 

2 

The  mean  square  refractive  index  fluctuation  < n > can  be  eliminated  in  favor  of  the 

2 

experimentally  measurable  Cn  that  appears  in  the  Kolmogorov  "two-thirds”  law, 

V »>  = Cn  °2/3  • 

Indeed,  using  the  power  series  expansion  of  the  Bessel  function  and  retaining 

only  the  first  two  terms  for  p/L ^ « 1 , the  Karman  structure  function  becomes: 

v2/3 


Dn(  P)  - 2 < n2>  (3) 

2L0 


Upon  comparison  with  the  Kolmogorov  equation, 


and  therefore: 


4=  -31  k7/6  l11'6 

in  perfect  agreement  with  the  result  of  an  earlier  analysis  performed  by  Tatarski. 
We  shall  now  make  use  of  the  asymptotic  behavior  of  the  Kummer  function, 


(1) 
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$ (a|c|x) 


with 


F(c)  . i c tt#  » a V (a)n(a-c^-l)n  , _.x-n 


r(c 


*h  <eie7*>*  l 


n! 


(-x) 


n=0 


N 


n=0 


(a>0  = i ; 


_ T( a+n ) ± 

(a)n  ’ 


C=1  if  Ira  x > 0 , € = -1  if  Irn  x < 0 ; 


- it  < arg  x < tt 


as  well  as  of  its  behavior  near  the  origin  to  extract  from  equation  D-45  information 
concerning  the  normalized  correlation  function  bB(^)  . One  can  thus  show,  after 

some  lengthy  algebra  that  as  long  as  0 does  not  exceed  tan  (/L/  X ) , 
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.5/6 

bA(  p)  - 1 - H/6  sin  9J 


for 


P « 


t an  8 

cos  0 


n/6 

25/3r(ll/6)  sin  (11/12  tt) 


( p cos 


5/3 

0) 


for 


t an  8 
cos  0 


« p « 


/XL 
cos  0 


(D-46) 


b »(?> 


j7/3-i2-a./6) (tl)7/6(,c„s  n,-7/3 

lire  11/6 )r(- 11/6)  Sin  (11/12  tt)  ' 


for 


P 


» 


/ XL 
cos  0 


It  follows  from  equation  D-46  that  the  normalized  correlation  function  bft(  p) 

starts  off  at  1,  decreases  as  one  minus  the  five -thirds  power  of  the  normalized  distance 
p cos  0 //XL  until  it  crosses  the  abscisa  somewhere  around  /xL/cos  0 , and  then 
approaches  zero  asymptotically  from  below  • 

An  exact  representation  of  the  correlation  function  b^(  p)  is  given  in  figure  4,  for 
0 = 10°  . 

The  correlation  length  of  the  logorithmic  amplitude  distribution  is  thus  of  order 
/XL 

^corr  ~ cos  0 
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For  0=0,  this  result  is  in  perfect  agreement  with  the  Tatarski  calculation  in  the 
/XL  « Lq  regime,  and  is  therefore  the  natural  generalization  thereof  away  from  trans- 
versal correlation.  As  0 nears  tt/2  the  result  of  equation  D-45  is  no  longer  accurate. 
In  fact,  for  9 = 9q  the  solution  of 

4tt  L0  = /XL 
X sin  0 cos  9 ’ 


2 2 

the  correlation  length  is  of  order  4n  Lq/  X sin  9 and  hence  values  of  p > 4tt  Lq/  X sin  9 

become  relevant.  It  can  be  argued,  however,  that  for  values  of  9 ^ 0q  the  correlation 
length  becomes 


4n  L0 

pcorr  ~ x sin  9 


(D-47) 


Indeed,  introducing 


H s SinJ 
4tt  l0 


we  can  still  use  equation  D-45  to  describe  the  behavior  of  the  normalized  correlation  func- 
tion for  ja  « 1 . When  |i  » 1 f we  can  return  to  equation  D-41  and  use  the  asymptotic 
behavior  of  the  Tricomi  function, 
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V n (a)n  (a_c+1\ 

Y (a|c|x)  - £ HI ' X 
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to  have : 
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2 < arg  x < tt 
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THETR  = 10,0 


• • • # 


FIG.  4:  THE  NORMALIZED  AMPLITUDE  CORRELATION  FUNCTION 


Corre  spondingly , 


K(z)  - - Ei  (-P2/4z) 

providing  for  the  correlation  function  a highly  oscillatory  behavior  in  the  |i  » 1 region. 

The  change  of  regime  manifested  by  5 ( p)  as  ja  goes  through  1,  indicates  that  the  Corre- 
ia 

lation  length  corresponds  to  values  of  u near  1 and  hence  that  equation  D-47  holds.  To 
obtain  the  correct  shape  of  the  correlation  function  b^(  p)  as  p passes  1 would  require 
summing  up  the  series  in  equation  D-41  without  any  approximation  on  the  Tricomi  functions 
2 

Y (n+l  jn+1/6  |z/Lq  ) . We  shall  restrain  from  attempting  to  do  so  here. 

For  9 = tt  /2  our  result  corresponds  exactly  to  the  one  obtained  in  1961  by  Chernov. 

5.  COMPARISON  WITH  EXPERIMENTAL  DATA 


The  value  of  the  variance, 


2 = . 31  C*  k7/6  L11/6  ; JTC  « LQ  , 


and  the  correlation  length 

/nr 


0 

Kcorr 


cos  9 5 0 5 < 90 


j / \L  « Lr 


4ttL0 
^ X sin  9 


9Q  < 9 ^ tt/2 


for  the  logarithmic  amplitude  fluctuation  field  of  a monocromatic  sound  plane  wave  propagating 
in  a turbulent  ocean,  are  our  main  result.  It  is  to  be  observed  that  they  depend  explicitly 
on  the  wave  number  k of  the  propagating  sound  and  on  the  distance  L between  the  source 

2 

and  the  receiver,  and  implicitly  through  Cn  on  the  environmental  conditions  prevailing 
in  the  region  of  propagation. 

We  shall  conclude  this  analysis  by  providing  some  information  concerning  the  compari- 
son of  our  result  at  9=0  with  adequately  taken  experimental  data . The  experiments  we 
have  in  mind,  were  conducted  in  1961,  and  then  again  in  1962,  1964,  from  two  scientific 
research  vessels  belonging  to  the  Academy  of  Sciences  of  the  USSR,  the  Sergei  Vavilov  and 
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(3) 

the  Peter  Lebedeev.'  ' Periodically  repeated  pulsed  signals  with  a tone  frequency  carrier 
f were  recorded  at  a definite  distance  from  the  transmitter.  The  depth  H of  the  trans- 
mitter and  receiver,  the  distance  L between  them,  the  pulse  duration  A t , and  their 
period  of  repetition  were  chosen  so  as  to  permit  time  separation  of  the  direct  and 
surface -reflected  signals.  As  a rule,  simultaneously  with  the  recording  of  the  acoustic 
signals,  the  fluctuations  of  the  velocity  of  sound  were  measured  at  the  depth  of  the  receiver. 
The  receiver  was  designed  to  operate  in  conjunction  with  a device  permitting  the  fluctua- 
tions of  the  refractive  index  to  be  measured  to  10"6  by  the  microphasometric  method.  This 
simultaneous  recording  procedure  provided  the  opportunity  of  analyzing  the  results  to  in- 
clude the  variation  of  the  statistical  characteristics  of  the  refractive  index  fluctuations  in 
the  course  of  the  experiment.  The  statistical  processing  of  the  recorded  fluctuations  of  the 
index  of  refraction  and  fluctuations  of  the  acoustic  pressure  level  of  the  direct  pulses  con- 
sisted in  calculating  the  time  structure  functions  D(  t)  according  to, 


D(  t) 


■ N-m 

I (zi+m 


i=l 


where  z^  is  the  value  of  the  ith  sampling  of  the  refractive  index  or  logarithmic  amplitude, 

N is  the  total  numbers  of  samplings  in  a recorded  segment  of  duration  T , At  = T/N  - 1 
is  the  partition  subinterval,  and  m = 1,  2,  ...  is  the  number  of  partition  subintervals 
At  in  the  interval  t . The  partition  subinterval  At  for  the  recordings  of  the  acoustic 
pulses  was  equal  to  their  period  of  repetition  AT  ; for  the  refraction  fluctuation  record- 
ings it  was  equal  to  0.05  second.  The  number  N was  between  the  limits  600  and  1,  200. 


The  transition  from  the  time  scale  t to  the  space  scale  o was  performed  accord- 
ing to  the  formula, 


0 = VT  . 

The  structure  functions  of  the  refractive  index  fluctuations  D^Cvt)  and  the  amplitude 
fluctuations  D^(vt)  were  used  to  ascertain  the  correlation  radii  and  mean-square  fluc- 
tuations of  the  refractive  index  and  of  the  logorithmic  amplitude.  The  principal  character- 
istics of  the  experiments  and  the  results  of  their  processing  are  shown  in  table  1 . 

2 

Figure  5 shows  the  normalized  structure  functions  Dn(  p)/2<n  > for  fluctuations  of 

the  refractive  index  for  the  recordings  obtained  in  the  1962  tests.  It  indicates  that  for 
values  of  0 < Lq  these  functions  on  the  average  fairly  well  satisfy  the  "two-thirds"  law  of 

Komogorov  and  Obukhov.  In  the  1964  experiments  direct  measurement  of  the  inhomogene- 
ities were  not  made,  but  there  existed  good  evidence  that  the  "two-thirds"  law  was  valid 
up  to  scales  of  at  least  60  miles. 
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TABLE  1:  EXPERIMENTAL  CHARACTERISTICS  AND  RESULTS 


TIME  AND  PLACE 
OF 
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• E 

200 

3.5 

9.5 

.36 

13.7 

.63 

1962  May,  Northwest  Atlantic 

40 

25 

1 

.4 

665 

6.3 

26.4 

.46 

25.4 

.28 

900 

7.4 

18.7 

.11 

29.0 

.11 

1150 

8.3 

14.1 

1.16 

4.7 

1.94 

1962  June,  Northeast  Atlantic 

35 

25 

1 

.4 
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7.5 

7.1 

1.58 

2.7 

1.98 
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6.0 

10.0 

2.22 

5.4 

2.66 
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5.4 

15.8 

.80 

11.0 

.80 

1962  July,  Sea  of  Norway 
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.4 
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FIG.  5:  THE  TEMPERATURE  STRUCTURE  FUNCTION 


2 

The  normalized  structure  functions  D^(  p)/2o^  obtained  from  processing  of  the 

recordings  of  a sequence  of  pulses  at  various  distances  are  constructed  as  a function  of 
p//xE  in  figure  6.  Notice  that  the  curves  for  the  structure  functions  roughly  coincide 
at  all  distances  and  intersect  the  saturation  level  between  values  of  the  normalized  dis- 
tance from  .9  to  1.3.  This  indicates  that  the  correlation  distance  is  of  order  /xE 
when  /xL  •«  Lq 


FIG.  6:  THE  AMPLITUDE  STRUCTURE  FUNCTION 


r 2 2 7 /6T% 

Finally  the  dependence  of  o^/.31  Cn  K ' J2  on  L is  depicted  in  figure  7.  The 
line  drawn  through  the  experimental  points  by  the  method  of  the  least  squares  almost  per- 
fectly coincides  with  the  theoretical  dependence  of 


The  foregoing  discussion  concerning  the  results  of  experiments  on  the  fluctuations  of 
the  sound  level  in  an  ocean  medium  containing  inhomogeneities  of  the  velocity  of  sound 
shows  that  a random  field  of  the  refactive  index  may  be  described  to  a satisfactory  approxi- 
mation by  means  of  structure  functions  satisfying  the  Komogorov -Obukhov  "two-thirds"  law 
with  external  inhomogeneity  scales  Lq  ranging  from  several  tens  of  centimeters  to  several 

tens  of  meters,  and  the  experimental  dependences  of  the  mean-square  and  radius  of  corre- 
lation of  the  acoustic  pressure  level  fluctuations  on  the  distance  are  in  good  agreement  with 
the  theoretical  predictions . 


100  200  500  1000  2000  5000  10000  lu(M) 

FIG.  7:  RANGE  DEPENDENCE  OF  CORRELATION  COEFFICIENT 
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